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Abstract 


We present a comprehensive study of the interaction of human serum albumin (HSA) with 
poly(acrylic acid) (PAA; number average degree of polymerization: 25) in aqueous solution. The 
interaction of HSA with PAA is studied in dilute solution as the function of the concentration of added 
salt (20 - 100 mM) and temperature (25 - 37°C). Isothermal titration calorimetry (ITC) is used to 
analyze the interaction and to determine the binding constant and related thermodynamic data. It is 
found that only one PAA chain is bound per HSA molecule. The free energy of binding AG^ increases 
with temperature significantly. AG^ decreases with increasing salt concentration and is dominated by 
entropic contributions due to the release of bound counterions. Coarse-grained Langevin computer 
simulations treating the counterions in an explicit manner are used study the process of binding in 
detail. These simulations demonstrate that the PAA chains are bound in the Sudlow II site of the HSA. 
Moreover, AG^ is calculated from the simulations and found to be in very good agreement with the 
measured data. The simulations demonstrate clearly that the driving force of binding is the release of 
counterions in full agreement with the ITC-data. 


Introduction 

Human Serum Albumin (HSA) is an important 
transport protein that interacts with substrates as 
e.g. fatty acids^ and pharmaceuticals in a spe¬ 
cific mannerP Also, transport by HSA plays an 
important role for renal clearance and the interac¬ 
tion of HSA with typical uremic toxins^^ as e.g. 
phenylacetic acid, indoxyl sulfate, and p-cresyl 
sulfate needs to be understood in detail.l^ Ure¬ 
mic toxins play an essential role in the excessive 
cardiovascular mortality and all cause mortality 
in patients with end stage renal disease. Elimi¬ 
nation of these substances is essential to reduce 
the high cardiovascular disease risk.l^ The ob¬ 
vious biological importance of this transport has 
led to a great number of studies of the interac¬ 
tion of HSA with various substrates, most no¬ 
tably by crystallography,!^ electron paramagnetic 
resonance spectro scopy ,!^ and isothermal titration 
calorimetry (ITC).!^^!^Thus, it is well established 
that the interaction of various substrates with HSA 
takes place mainly on the Sudlow I and the Sud¬ 
low II site.!^ ITC has been used to elucidate the 
thermodynamics of this interaction and by this the 
number of subst rates adsorbed onto the surface of 
the protein.!^^^i®i3 

Charge-charge interaction plays an essential role 
in the process of adsorption!iSH22l and a number 
of ITC-studies explore the dependence of the ad¬ 
sorption constant on ionic stren gth.!^^^!^ Recently, 
Jankowski and coworker s!25l2il have demonstrated 
that raising the ionic strength in the infusion fluid 


leads to an improved clearance of protein-bound 
toxins (PBT). Thus, raising the concentration of 
NaCl to 600 mM led to a significantly better re¬ 
moval of of uremic toxins as e.g. phenylacetic 
acid. This result points clearly to the central im¬ 
portance of Coulombic interactions for the bind¬ 
ing strength of such toxins to HSA and to proteins 
in general. Charged toxins as phenylacetic acid 
are difficult to remove by conventional dialysis and 
exploring their interaction is of central importance 
for devising improved techniques for dialysis.!^ 

Here we analyze the role of the interaction of 
charged substrates with HSA as the function of 
ionic strength and temperature. As a substrate we 
chose a short polyelectrolyte, namely poly(acrylic 
acid) (PAA) with 25 repeating units only. This 
substrate allows us to explore the Coulombic in¬ 
teraction of HSA with charged molecules in gen¬ 
eral. In addition to this, synergetic effects of ad¬ 
jacent carboxyl groups onto the binding constant 
can be elucidated. At the same time, PAA can 
be regarded as a model of charged toxins with 
molecular weight above 500 g/Mol.!^] ITC is used 
to obtain the binding constant and the number of 
bound PAA-molecules per protein and the temper¬ 
ature is varied between 25°C and 37°C. A first 
study of the interaction of polyelectrolytes with 
proteins by ITC has been presented by Schaaf et 
al. who demonstrated the general suitability of this 
method for the study of protein-polyelectrolyte in¬ 
teraction.!^ 

The interaction of long polyelectrolytes with 
proteins in aqueous solution has been the subject 
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of longstanding researc h that has led to an enor¬ 
mous literatureJ ^^I^^ I ^^ I Proteins can form com¬ 
plex coacervates with polyelectrolyte of opposite 
charge in aqueous solution and the strength of in¬ 
teraction is mediated by the ionic strength in the 
system® If the ionic strength of the system is low 
enough, interaction may take place even on the 
“wrong side“ of the isoelectric point, that is, pro¬ 
teins associate with poly electrolytes of like charge. 
In many cases the formation of complexes between 
the protein and the polyelectrolyte is followed by 
precipitation and phase separ ation that may also 
lead to non-equilibrium states.®^ Here we use a 
very short polyelectrolyte and low concentrations 
to avoid multiple interactions and phase separa¬ 
tion. The results from the present experiments, 
however, may be directly utilized for a better un¬ 
derstanding of complex coacervates. 

Evidently, protein crystallography cannot be 
used to elucidate the structure of the complex be¬ 
tween HSA and PAA. First of all, even in the 
bound state the non-bound segments of the poly¬ 
electrolyte will explore their conformational free¬ 
dom in solution and the entropic contribution de¬ 
riving therefrom will be an important part of the 
resulting change of free energy. Moreover, crystals 
of HSA will certainly not accommodate substrates 
of the size of a poly electrolyte having 25 units. In 
order to make progress on a detailed structural in¬ 
vestigation, we employ coarse-grained Langevin 
Dynamics (LD) computer simulations with an im¬ 
plicit solvent whereas the co- and counterions are 
treated in an explicit manner.®® protein is 
treated within the - Go model, that is, each 
amino acid of HSA is modelled by a single bead 
bearing a charge or not. The polyelectrolyte PAA 
is also treated as a coarse-grained charged poly¬ 
mer. The combination of these models allows us 
to elucidate the details of the complex of PAA and 
HSA in solution with full structural resolution. In 
addition to this, LD-simulations allow us to calcu¬ 
late the free enthalpy of binding AGt, and to com¬ 
pare these data with measured values. Thus, the 
combination of ITC with computer simulation can 
be used to elucidate structural and thermodynamic 
details that are available by no other method. 


Experimental 

Materials. Polyacrylic acid (PAA) with 
Miv=1800 g/mol was purchased from Sigma- 
Aldrich (Schnelldorf, Germany) and used after 
several weeks of dialysis to match pH without 
changing ionic strength in the system. Human 
Serum Albumin (HSA) was also purchased from 
Sigma-Aldrich (lyophilized powder. Fatty acid 
free. Globulin free, 99%) with molecular weight 
calculated Mvt?=66400 g/mol and its purity verified 
by sds-gelelectrophoresis. The buffer morpholin- 
N-oxide (MOPS) was purchased from Sigma- 
Aldrich and used as received. 

Isothermal Titration Calorimetry 

ITC experiments were performed using a VP- 
ITC instrument (Microcal, Northampton, MA). 
All samples were prepared in a pH 7.2 buffer so¬ 
lution using 10 mM MOPS and 10 mM, 40 mM, 
60 mM and 90 mM NaCl to adjust ionic strength. 
All samples were dialyzed against buffer solu¬ 
tion with according pH and degassed prior to 
experiment. For dialysis, the dialysis-system 
Float-a-Lyzer by Spectrum Labs with molecular 
weight cut-off (MWCO) 500-1000 Da for PAA 
and MWCO 20 kDa for HSA were used respec¬ 
tively. The samples were thermostatted and the in¬ 
strument stabilized for 1 h to ensure thermal equi¬ 
librium and stability of the system. A total of 
298 jiL PAA solution was titrated with 75 succes¬ 
sive 4 jih injections, stirring at 307 rpm and a time 
interval of 300-350 s between each injection into 
the cell containing 1.4 mL protein solution. The 
concentration of PAA and HSA were 1.3 g/L and 
0.9 g/L respectively. We choose these low concen¬ 
trations in order to study the interaction of single 
PAA-chains with only one HSA-molecule. More¬ 
over, possible complications by coacervate forma¬ 
tion are circumvented in this way. The experi¬ 
ments were performed at 25, 27.5, 30, 33, 37 °C 
and ionic strengths 20 mM, 50 mM, 70 mM and 
100 mM. 

As a first step of the ITC data analysis, the inte¬ 
gration of the measured heat Q over time is carried 
out to obtain the incremental heat AQ as a function 
of the molar ratio x between polyelectrolyte and 
protein. For each experiment, the heat of dilution 
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of PAA were measured separately by titrating the 
according polyelectrolyte into blank buffer solu¬ 
tion and subtracted from adsorption heats. After 
correction, the resulting binding isotherm are fitted 
using a supplied module for Origin 7.0 (Microcal). 


Data analysis. 

For the analysis, we chose the single set of in¬ 
dependent binding sites (SSIS) model to fit all 
data.I^This model is based on the Langmuir equa¬ 
tion, which assumes an equilibrium between the 
empty adsorption sites, the no. of proteins in solu¬ 
tion and the occupied adsorption sites. This leads 
to the binding constant Ki,: 


0 

(1 —&)CpAA 

where 0 denotes the fraction of sites occupied by 
the poly electrolyte and cpaa the concentration of 
free polyelectrolyte. Since only the total concen¬ 
tration of PAA in the solution is known, cpaa 

is connected to the as follows: 


Furthermore using the following two equations, 
the entropy A5/, can be either calculated for one 
temperature or extracted from its temperature de¬ 
pendence: 


AGb = AHjTC-TAS (5) 


dAGb 

dT 


= —^b 


( 6 ) 


Theoretical Methods 

Computer simulation model and parameters 

In our simulations we employ an implicit-water 
coarse-grained (CG) model, where each of the 
amino acids, PAA monomers, and salt ions is ex¬ 
plicitly represented by a single interaction bead. 
Hence, the water is modeled by a dielectric back¬ 
ground continuum while the salt is explicitly re¬ 
solved. The dynamics of each of the beads in our 
simulations is governed by Langevin’s equation of 
motion!^ 

r ■ H T' 

-|-V,C/+ (7) 


^PAA — ^PAA T AfQCpfot (2) 

where N is the number of free binding sites and 
Cproi the total protein concentration in solution. 
Following these equations, the binding number 
N, binding affinity Kj, and the overall enthalpy 
change measured AHjpc can be obtained by fitting 
the isotherm. In the following, we will show that 
only one PAA is adsorbed onto HSA in all cases. 
Therefore the parameter N was fixed to unity for 
all subsequent fits. In this case the interaction of 
PAA and HSA can be formulated as a conventional 
chemical equilibrium: 


CPAa{(^HSA ~ <^PAA-HSA 

where denotes the total protein concentra¬ 
tion in the solution and cpaa-hsa the concentration 
of PAA-HSA complex. Furthermore, the binding 
free energy AGy can be derived by 

AGt = -RT ■ InKb (4) 


where m,- and are the mass and friction constant 
of the ith bead, respectively. U is the system po¬ 
tential energy and includes harmonic angular and 
bonded interactions between neighbouring beads 
in HSA and PAA, dihedral potential in the case of 
the HSA, and interatomic Lennard-Jones (LJ) be¬ 
tween all non-neighbouring beads, including ions. 
Coulomb interactions govern the electrostatic pair 
potential between all charged beads. The random 
force Ri{t) is a Gaussian noise process and satis¬ 
fies the fluctuation-dissipation theorem 

{Ri{t)-Rj{t')) =2mi^ikBT5{t-t')5ij. (8) 

The simulations are performed using the GRO- 
MACS 4.5.4 software package.^^A leap-frog al¬ 
gorithm with a time step of 2 fs is used to inte¬ 
grate the equations of motion. A cubic box with 
side lengths of L = 30 nm is employed and period¬ 
ically replicated to generate a quasi-infinite system 
in the canonical ensemble. The Langevin friction 
is =1.0 ps“^ that dissipates the energy at con¬ 
stant temperature T. Center of mass translation of 
the system is removed every 10 steps. The cutoff 
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radius is set to 3.0 nm to calculate the real-space 
interactions while Particle-Mesh-Ewaldl^ (PME) 
is implemented to account for long-range elec¬ 
trostatics. The PME method is computed in the 
reciprocal space with a EFT grid of ~ 0.12 nm 
spacing and a cubic interpolation of fourth-order. 
The solvent is modelled as a continuous medium 
with a static dielectric constant of = 73.4 and 
78.2 for temperatures T = 25°C and 37°C, respec¬ 
tively. All beads have mass m,- = 1 amu, diameter 
Ojj = 0.4 nm, energy well ejj = 0.1 keT and in¬ 
teger charges g = 0, -i-l or -1 e. The mass was cho¬ 
sen artificially low to enhance orientational fluctu¬ 
ations and sampling. Clearly, equilibrium proper¬ 
ties, as investigated in this work, are not affected 
by any reasonable mass choices as long as the sim¬ 
ulations are ergodic. 

The protein sequence for the HSA is provided 
by PDB databank: ID=lN5U.I^Every amino acid 
is described by a single bead positioned at its Ca 
atom. The native structure of the protein is main¬ 
tained by a Go-model like force field as provided 
from t he SM OG webtool for biomolecular simu¬ 
lations All beads corresponding to basic and 
acidic amino acids are assigned a charge of +le 
and -le, according to their titration state at physi¬ 
ological pH = 7.4, that is, ARG and LYS residues 
are assigned a charge +le, while ASP and GLU are 
-le, and HIS is neutral. With that the net charge of 
the simulated HSA is -14e. 

A single flexible PAA polyelectrolyte is mod¬ 
eled in a coarse-grained fashion as a sequence of 
Nmon freely jointed beads. Each bead represents 
a monomer with a radius Gu and carries an elec¬ 
tric charge of — le. The PAA monomers are con¬ 
nected by a harmonic bond potential with an equi¬ 
librium bond length ^mon = 0.4 nm and a force 
constant A:nion = 4100 kJ mol“^ nm“^. The flex¬ 
ibility of the PAA chain is defined via a harmonic 
angle potential in which the angle between a triplet 
of monomer is ^ = 120° and the force constant is 
k(i) = 418.4 kJ mol“^ rad“^. As in the experiments 
we consider short PAA chains with Amon = 25 
monomers. The simulation box contains one HSA, 
while the amount of PAA and ions are character¬ 
ized by the molar ratio x = cpaaIchsa, ranging 
from 1 to 10, and salt concentration Csait, ranging 
from 20 to 100 mM, respectively. 


Binding and free energy calculations 

The stoichiometry, in our case the average num¬ 
ber of bound PAA chains on one HSA, can be de¬ 
termined through a calculation of the normalized 
density distribution function g{r) = c{r)/c’pj^, 
where r is the distance between the centers-of- 
mass (COM) coordinates of the HSA and PAA and 
is the PAA bulk concentration. Integration 
of the g(r) further leads to the PAA coordination 
number 

n(r) = 4;rc^Xi / g{r')r'^dr'. (9) 

Jo 

which describes how many PAA chains are bound 
on average at a distance r. For each molar ratio 
we simulated 120 ns in order to calculate the equi¬ 
librium coordination (binding) numbers of PAA to 
HSA. To quantify the number of bound and re¬ 
leased ions upon complexation, we count the av¬ 
erage number of ions Nj, i = ±1, that are bound 
(’condensed’) on the PAA chain or on the positive 
protein patches, respectively, and make a compari¬ 
son before and after PAA/HSA association. An ion 
is defined as ’condensed’ if it is located in the first 
binding layer, that is closer than a cut-off distance 
rs = 0.5 nm from the charged bead, while double¬ 
counting in overlapping volumes is avoided. The 
average is over long (ca. 30 ns) trajectories in the 
fully separated and stable bound states. To calcu¬ 
late the potential mean force (PME) between the 
HSA and the PAA we employed steered Langevin 
Dynamics simulations using the pull code as pro¬ 
vided by GROM ACS.Here, the center of mass 
of the PAA is restrained in space by an external 
time-dependent force. This force is applied as 
a constraint, i.e. harmonic potential, and moved 
with a constant pulling velocity Vp to steer the par¬ 
ticle in the prescribed direction.!^ After several 
test runs, the pulling rate Vp = 0.1 nm/ns was cho¬ 
sen and a harmonic force constant K = 2500 kJ 
mol“^ nm“^. The simulations were performed for 
~100 ns. Given the pulling speed above, this sim¬ 
ulation time is required to bring the two macro¬ 
molecules from a separated state (r ~ 11 nm) to a 
final state (r ~ 1 nm). The standard deviation was 
calculated by standard block averages to specify 
the statistical error. The friction force F = —m^jVp 
was subtracted from the constraint force and the 
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result averaged within a specific interval of dis¬ 
crete spacing Ar to obtain the mean force of the 
interaction potential. According to the simulation 
setup, the mean force was integrated backwards to 
get the potential of mean force (PMF). Because 
the CPPM were radially constraint in 3D space, 
the PMF has to be corrected for translational en- 
tropyl^by 

G(r) = G^(r)-2kBTln(r), (10) 

where G^(r) is the integrated mean force. The 
binding affinity of the PAA can then be defined as 


the free energy value at the global minimum of the 
PMF in the stable complex as 

AG'*"” = G(rm,'„) — G(oo), (11) 

where the reference G{°°) is set to zero. How¬ 
ever, before making a comparison with the exper¬ 
iment, we had to consider that AG^^^ provided by 
the experiment is defined as a standard free en¬ 
ergy, which refers to the standard binding volume 
Vq = \/C^ of one liter per mol.l^ Hence, the bind¬ 
ing constant K}, generated from experiment is for- 




Figure 1: ITC data of adsorption of PAA upon HSA and the corresponding heats of dilution of PAA at pH=7.2, T=37°C 
and a) l=20 mM, b) l=50 mM. c) Binding isotherm corrected for the heat of dilution at 37°C and l=20 mM & 50 mM. 



Figure 2: Determination of the parameter N, i.e., the number of PAA molecules adsorbed on on HSA molecule is 
shown for one temperature 37°C at a) i=20 mM and b) 50 mM. Each iine shows a fit with different vaiues of fixed N 
marked by different colors. 
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mulated as In our simulations 

we average the accessible volume I 4 of the COM 
of the PA A in the bound state. 

As a result, the standard binding free energy 
from the simulation can be obtained as^^ 

AGf ” = , (12) 

with a term AG'^'^" = —ksT ln(C*^V),) is the entropy 
correction arising from the accessible volume of 
the COM of the PAA in the bound state. 

Results and Discussion 

We performed a systematic series of ITC exper¬ 
iments comprising four ionic strengths and five 
different temperatures ranging from room tem¬ 
perature (25°C) to the physiological temperature 
(37°C). The experiments were performed at pH 
7.2 in buffer solution. This pH is well above the 
isoelectric point of the protein HSA used for this 
experiment leading to a net effective charge of -14. 
PAA is a weak acid with a pK^ of 4.5 and hence 
a net charge of -25 at pH 7.2. Thus, we study the 
adsorption on the "wrong side" of the isoelectric 
point where both the protein as well as the poly¬ 
electrolyte is charged negatively. 

ITC is certainly the method of choice for the de¬ 
termination of th e adsorp tion constant of a given 
substrate to HS A.^^ElIS] jjowever, the concentra¬ 
tions of the protein and PAA are small and the 
evolved heat will be concomitantly small. Hence, 
all effects leading to spurious heat signals must be 
considered in detail and carefully excluded. The 
main problem is the adjustment of the same pH 
and ionic strength in both the solution of the pro¬ 
tein and of the polyelectrolyte. This is done by 
extensive dialysis which turned out to be decisive 
for obtaining meaningful ITC-data. Since the con¬ 
centrations of both PAA and HSA are low, all spu¬ 
rious effects leading to a heat signal must be con¬ 
sidered. Since PAA-solutions are added in small 
portions to the solution of HSA, the heat of dilu¬ 
tion of PAA must be determined carefully and sub¬ 
tracted from the raw signal. Fig. [T] and Fig. SI of 
the SI displays the raw ITC-signals of PAA onto 
HSA (black curves and points) as well as the heat 
of dilution of PAA (green points). The signal is 


weakly endothermic at 1=20 mM but exothermic 
at 1=50 mM (see Fig.[T]). Dilutions are in all cases 
exothermic and the effect becomes stronger with 
increasing salt content as expected (blue curves 
and points). For higher ionic strength, the heat of 
dilution has a dominant effect on the overall signal 
and determines the sign of the signal at I>20 mM. 
For data analysis, the heats of dilution are sub¬ 
tracted from the heats of adsorption prior to fitting. 
Here, special attention must be paid to this step, as 
for some cases there remains a constant residue af¬ 
ter subtraction of the heat of dilution. Even though 
this offset signal is very small and usually less 
than 0.1 kcal/mol, it cannot be neglected due to the 
small overall heat. We assign this offset to a slight 
mismatch of the pH or salt titrant and the solution 
in the cell. In order to take this effect into account, 
a flat background was fitted to all isotherms after 
the first step of subtraction and used to correct the 
data in the second step. 

The panel on the right-hand side of Fig. dis¬ 
plays a set of typical results thus obtained. Evi¬ 
dently, the heat of adsorption is positive and we 
find this for all conditions under consideration 
here. Hence, the driving force for the process of 
adsorption must be entropic. This point will be 
discussed in more detail below and is well borne 
out from the simulations, too. 

In order to obtain the number N of PAA- 
molecules bound to one HSA-molecule, the data 
were first fitted using the SSIS model as described 
in Section 2.1.1. Fig. shows a comparison of the 
parameter N for the two data sets as in Fig. The 
colored curves showing different fixed values of N 
reveal that the data clearly justify A = 1 as the best 
choice for fitting. This observation is true for any 
other sets of data. Deviations from N = 1 are not 
significant and we can safely assume A = 1 in all 
subsequent analysis (see eq. ?? in Section 2.1.1). 

Strength of interaction as the function 
of temperature 

Fig. presents two series of temperature depen¬ 
dency studies for the ionic strengths 1=20 mM and 
50 mM. For better clarity, only 3 temperatures 
are displayed in these graphs. The data for two 
more temperatures are displayed in Figure SI of 
the Supporting Information. The data taken at both 
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Figure 3: Effect of temperature. The integrated heats Q for adsorption of PAA on HSA at temperatures between 25°C 
and 37°C for a) 1=20 mM and b) 1=50 mM and the according fits are displayed. For better clarity, only 3 temperatures 
are displayed, c) Van’t Hoff analysis of the dependence of the adsorption constant on temperature. Data points are 
derived from the fits of the ITC-data shown in a) and b). The crosses correspond to 1=20 mM and the empty squares 
to 1=50 mM. 


ionic strength reveal a significant increase of en¬ 
thalpy with increasing temperature from 25°C to 
37°C. This effect is more pronounced for 1=20 mM 
than at 1=50 mM. Additionally, the overall en¬ 
thalpy of adsorption becomes weaker with increas¬ 
ing salt content, which points directly to the im¬ 
portance of electrostatic interaction on the bind¬ 
ing of PAA onto HSA. All data are very well de¬ 
scribed by the model assuming N = 1. Data taken 
at higher salinity are more noisy but the raise of 
the signal with temperature is clearly discernible. 
Since there is no plateau in the ITC signal, the pa¬ 
rameter AHjtc might be overestimated by the fits. 

The results of these fits are listed in Tabled with 
AHjtc and Kj, as fitting parameters, N being fixed 
to unity (see above). The free energy of binding 
AGh was calculated from the fitting parameter 
using equation ??. 

The temperature dependence of the binding can 
now be analyzed according to van’t Hoff’s law: 

(SS),-? ™ 

Fig. c) shows the resulting van’t Hoff plot. A 
linear correlation between logarithm of the bind¬ 
ing constant Kj, and the inverse temperature is seen 
within the limits of errors. The binding enthalpy 
Mil, can be obtained from the slope of the lin¬ 
ear fit and ASj, from the intercept. The resulting 
data are gathered in Table [Tj In general, the values 


of AsRjtc larger than the data resulting from 
the van’t Hoff analysis. Similar finding have been 
made in a recent study of the interaction of pro¬ 
teins with charged microgels.^^ Reasons may be 
sought in additional processes as e.g. the hydra¬ 
tion of freed counterions that are not directly cou¬ 
pled to the process of binding (see below). Also, 
the heat of adsorption taken directly from the ITC 
data might be slightly overestimated (see above). 

Dependence on ionic strength. 

To study the dependence of the binding process 
on ionic strength, we carried out two more exper¬ 
iments at 37°C and 1=70 mM and 100 mM (see 
Fig.|^ raw data are shown in Figure S2). With in¬ 
creasing ionic strength, the measured enthalpy ap¬ 
proaches zero and the ITC method reaches its in¬ 
strumental limits. As described above, we fix the 
parameter N for both salt concentrations to unity. 
Table [^gathers all data obtained from these exper¬ 
iments. 

The data exhibit a very consistent decrease of 
binding constant Kj, with increasing salt concen¬ 
tration (see Tab. [T]). This observation combined 
with the fact that only about one PAA molecule is 
adsorbed on the HSA leads to the conclusion that 
the driving force of the interaction is an attractive 
electrostatic potential between the negative PAA 
with patches of positive charge on the surface of 
HSA molecule. These patches are known to act 
as multivalent counterions for the poly electrolyte. 
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Table 1: Overview of thermodynamic parameters for all fitted isotherms for the temperature series between 25°C- 
37°C and ionic strengths between I = 20 mM -100 mM. As discussed in Sec. 3, all data were fitted with fixed N =1. 

ASh and AHh were calculated according to eq. ??,?? and ?? respectively. Entropys for 70 mM and 100 mM 
were calculated using eq. ?? at 37°C. 


Ionic Strength (mM) T (°C) 

AHjtc 

(kJ/mol) 

Kb • 10"^ 
(mol“^) 

(kJ/mol) 

ASb 

(kJ/mol/K) 

AHb 

(kJ/mol) 


25 

16.4±0.3 

7.7±0.5 

-27.9±1.3 




27.5 

32.8±0.6 

5.8±0.3 

-27.4±1.4 



20 

30 

34.0±0.5 

6.9±0.3 

-28.1±1.1 

0.17±0.01 

15±4 


33 

45.6±0.9 

8.1±0.5 

-28.8±1.6 




37 

53.4±1.0 

8.3±0.5 

-29.2±1.5 




25 

13.6±1.1 

1.3±0.2 

-23.4±0.5 




27.5 

14.4±0.9 

1.6±0.2 

-24.2±0.3 



50 

30 

24.8±0.7 

1.7±0.1 

-24.6±0.1 

0.27±0.02 

44±8 


33 

27.6±0.9 

1.8±0.2 

-24.9±0.2 




37 

26.4±0.6 

2.6±0.1 

-26.2±0.1 



70 

37 

21.9±0.7 

l.OiO.l 

-23.6±0.5 

0.12 (37°C) 


100 

37 

35.5±11.6 

0.3±0.1 

-20.7±1.2 

0.10 (37°C) 



Binding of such a polyelectrolyte to such a patch 
thus leads to a release of its counterions.SSfiZljhis 
“counterion release force“ was considered many 
years ago by Record and Lohman who predicted a 
linear correlation between the logarithms of bind¬ 
ing constant and salt concentration^^ 


dlnKi, 

dlncsait 


(14) 


where K/, is the binding constant, Cgait the salt con¬ 
centration and Hion the number of counterion re¬ 
lease upon adsorption. This behavior is observed 
for the present data as well if the ionic strengths is 
above 20 mM (see Fig. |^. Application of eq. ?? 
to the data in Fig. [^yields Anion ~ 2.9 ±0.5, that 
is, approximately 3 ions are released upon bind¬ 
ing of one PAA-molecule to a HSA molecule. The 
deviation from linearity for low ionic strength has 
been observed as well by Dubin et al. for the in¬ 
teraction between Bovine Serum Albumin and the 
poly anion Heparin at pH 6.8.1^ Here, electrostatic 
interactions become long-ranged and the relative 
contributions of the counterion release mechanism 


significantly decrease. 

Extrapolation of the linear fit reveals that inter¬ 
action still persists at physiological ionic strengths 
and temperature. At 1=150 mM and 37°C, we 
derive from our plot a finite binding free en¬ 
ergy of around -17 kJ/mol and it only decays to 
small values at around 750 mM. This concen¬ 
tration has been found to be necessary in dialy¬ 
sis to remove protein-bound uremic retention so¬ 
lutes.!^ Evidently, there is still some binding un¬ 
der physiological conditions and much higher salt 
concentrations are needed to remove the toxins 
from the surface of HSA. These results show also 
that ITC-experiments can be deceptive when en¬ 
thalpies of different reaction in the system com¬ 
pensate in such a way that they vanish. Thus, in 
the present case interaction still exists under phys¬ 
iological conditions but does not lead to measur¬ 
able enthalpies. 
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Figure 4: Effect of ionic strength. Isotherms are shown 
for a series of ionic strengths ranging from 1=20 mM- 
100 mM at 37°C. All fits have been done with N — 
The thermodynamic data resulting from these fits are 
listed in Table[^ 



Figure 5: Linear dependence of the binding constant 
versus ionic strength as listed in Tablejilare depicted in 
logarithmic scales according to eq.??.l^ 

Comparison with computer sim¬ 
ulations 

All simulations presented in the following are 
using an implicit-water coarse-grained (CG) 
structure-based model, where the CG protein 
is held in its native structure by semi-empricial 
force-fields.l^ Here the amino acids, the PAA 
monomers, and the salt ions are modelled explic¬ 
itly on a single bead level. Water is modelled by 
a dielectric background continuum while the salt 
ions are explicitly resolved. Similar methods have 
been used repeatedly to stu^ e.g., protein fold- 
ingl^ and the pair potentially between proteins. 


They provide a reasonable picture of the interac¬ 
tions of a given molecule with the amino acids 
localized at the surface of the protein because they 
retain native structure, thermal motion, and ion- 
induced mechanisms in a well-resolved fashion. 
Hence, this type of simulations provides a full mi¬ 
croscopic picture of the interaction of PAA with 
HSA, in particular the identification of the sites 
where PAA docks on. Moreover, the simulations 
can be used to obtain realistic free enthalpies of 
binding as will be further shown below. Our com- 



Figure 6: (a) Representative computer simulation snap¬ 
shot of the total HSA-PAA complex, (b) Magnification 
of the binding site: the PAA (yellow string of beads) is 
bound near the Sudlow II site. The amino acid beads 
that directly participate in the binding (defined by be¬ 
ing within 0.5 nm distance to the PAA on average) are 
depicted by the opaque spheres. The rest of the HSA 
structure is distinguished by a transparent surface plot. 
Electrostatically neutral HSA beads are colored white, 
positive beads are green, and negative beads are red. 

puter simulations demonstrate that the HSA binds 
only one PAA chain, independent of temperature, 
salt concentration, and molar ratio in the consid¬ 
ered parameter ranges. Hence, it reconfirms the re¬ 
sult obtained previously by ITC. A representative 
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simulation snapshot of the bound complex with 
one PAA is presented in Fig. From a thorough 
screening of our simulation trajectories it emerges 
that this structure is highly reproducible and as¬ 
sumed in 80% of the simulation time. It is hence 
a highly stable and probable configuration. Ad¬ 
ditional analysis reveals that the PAA chain spans 
the sub-domains II A, III A, and III B, involv¬ 
ing the Sudlow II binding site. As expected for 
a negatively charged polyelectrolyte we find that 
it favorably binds positively charged amino acids, 
arginine (R) and lysine (K) at positions R410, 
R484, R485, R413, R538, K541, K199, K195, see 
also the green opaque spheres in Fig. This is 
certainly a central result of the present analysis 
inasmuch as it defines the precise location of the 
binding of a highly charged molecule as PAA. 



Figure 7: Running coordination number n{r) of the PAA 
chains around the HSA versus their centers-of-mass 
distance r at a temperature of 25° C, salt concentra¬ 
tion of 20 mM, and molar ratios x = cpaa/chsa = 1,2, 
6, and 10, see legend. In all cases roughly one PAA 
chain (horizontal black dotted line) is bound to HSA. 

To demonstrate that only one PAA forms a com¬ 
plex with HSA, the running coordination number 
n{r) of the PAA around HSA is shown in Fig. 
for molar ratios x = 1, 2, 6, and 10. The quan¬ 
tity n{r) is the total number of PAA-molecules 
around a given HSA as the function of distance. 
The plateau of the curves between separations of 
r = 4 and 6 nm just above the average value of 
the HSA radius is a proof that the binding num¬ 
ber does not exceed one, irrespective of the molar 
ratio. Only at larger distances n{r) is increasing 
beyond unity since here the entire solution is ex¬ 


plored. We find qualitatively similar results for the 
other investigated salt concentrations and temper¬ 
atures. This finding again is in direct agreement 
with the results of the ITC-experiments (see also 
the normalized density profiles in Fig. S5). 

It is furthermore interesting to see the temporal 
evolution of the complex of HSA with PAA. The 
PAA chain slides along the Sudlow II site much 
in a way of a threading through an orifice. A se¬ 
ries of snapshots combined with a movie can be 
found the in Fig. S4 of the supplementary infor¬ 
mation. At the one hand, this fact demonstrates 
the strong binding of PAA by this side. On the 
other hand, the threading through this site leads to 
a strongly increased number of configurations of 
the complex and thus increases the entropy of the 
complex. This fact certainly leads to the binding 
of the PAA-chain at the Sudlow II site and not on 
other positive patches on the surface of HSA. 

Our simulations allow us to calculate the free 
energy profiles (potential of a mean force) along 
the HSA-PAA distance coordinate. Examples for 
this interaction free energy G(r) between a sin¬ 
gle uncomplexed HSA and one PAA at two salt 
concentrations is presented in Fig. For larger 
distances of approach, r ~ 7 nm, a small repul¬ 
sive barrier can be observed stemming from the 
monopole charge repulsion which dominates for 
large separations as expected. The barrier de¬ 
creases and shifts slightly to shorter distances with 
higher salt concentration. At about r ~ 6 nm 
the onset of a strong attraction takes place until 
a global minimum is observed at closer approach 
at about r = 2 nm. The onset of attraction hap¬ 
pens right at values comparable to half of the con¬ 
tour length of the PAA chain at which one of its 
ends is first able to get in contact with the HSA 
surface. We never found a stable free energy min¬ 
imum for the adsorption of a second PAA. This is 
due to a too strong monopole charge repulsion and 
the covering of the high-potential binding spot by 
the firstly bound PAA. 

For the stable HSA/PAA complex, the binding 
free energy AG''"” can be calculated from the dif¬ 
ference of the global minimum and the reference 
free energy at large distances (horizontal lines in 
Fig. [^. The values of the simulation binding free 
energy, corrected to yield the standard free energy 
of binding from the simulations AG^™ (see Meth- 
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ods), are summarized in Table for various salt 
concentrations and temperatures. We find a very 
good agreement for all systems with largest de¬ 
viation of only 13% for the case of 50 mM salt 
at 37° C temperature. As in the experiments, the 
binding affinity decreases in the simulations with 
higher salt concentrations and increases with ris¬ 
ing temperature. The highly quantitative descrip¬ 
tion by the simulations, however, is actually some¬ 
what surprising given the simplicity of the under¬ 
lying model and the neglect of hydration effects 
and should be discussed with caution. However, 
we take it as a strong indication that relatively 
generic electrostatic interactions rule the complex- 
ation process and hydration contributions (such as 
hydrophobic or van der Waals (vdW) attractions) 
are rather small. 



Figure 8: Free energy profile (or potential of mean 
force) G(r) between the PAA and the HSA versus their 
centers-of-mass distance r at a temperature of 25° C 
and for 20 mM (red) and 50 mM (blue) salt concentra¬ 
tions. The binding free energy AGgim derived from the 
simulation can be read off as the difference between the 
zero free energy reference state at far separation (hori¬ 
zontal black dotted line) and the global minimum repre¬ 
senting the bound state (horizontal blue and red dotted 
lines). 

In order to test our hypothesis introduced above 
that the binding free energy is essentially domi¬ 
nated by the entropy of released counterions from 
the PAA chain and/or a positive patches on the 
HSA, we have counted the number of released ions 
upon complexation. In brief, we define an ion as 
’condensed’ if it is located in the first bound layer 
near a charged HSA or PAA monomer, defined by 


a cut-off radius of 0.5 nm. Hence, the number of 
released ions is calculated as the difference of the 
average number of condensed ions in the fully sep¬ 
arated and the stable bound states (see Fig. S5). 
For a temperature of 25° C and 20 mM and 50 mM 
salt concentrations, our analysis shows that on av¬ 
erage indeed 2.5 condensed ions are diluted away 
into the bulk upon complexation. This is indeed in 
good agreement with the Record-Lohman-analysis 
of the experiment data discussed above (cf. the 
discussion of Fig. |^. 

Deeper inspection shows that 2 of those ions 
come from the PAA, at which they where bound 
in a high density state (see Fig. S6 for ionic pro¬ 
files around PAA monomers). If we now just con¬ 
sider the PAA-condensed ions and their average 
concentration in the bound state ~ I.5± 

0.5 mol/1, this implies that a favorable entropy con¬ 
tributions of about ~ ksTln^Cdense/cs) = 4.3 ±0.4 
and 3.4 ±0.4 ksT is gained per ion upon its re¬ 
lease into 20 mM and 50 mM bulk concentrations, 
respectively. The total release free energies es¬ 
timated by this analysis are thus roughly -21 ±2 
and -17±2 kJ/mol for 20 mM and 50 mM salt, re¬ 
spectively, which is close to the binding free en¬ 
ergies from both experimental data and from sim¬ 
ulations. Hence, the binding of PAA is to a great 
part ruled by a counterion release mechanism and 
entropy. We note, however, that the matching of 
these numbers may be fortunate since other non- 
negligible interactions such as (repulsive) chain 
entropy, vdW attractions, and multipolar charge 
interactions beyond the bound ion layer (that is, 
from screening ions), all present in both simula¬ 
tion and experiment, have been neglected in this 
simple counterion release concept. The present 
comparison with experimental data, however, in¬ 
dicates that these contributions are of comparable 
magnitude and cancel each other roughly for the 
present system. 

Evidently, we obtained the leading contribution 
of the ions directly condensed on the PAA chain. 
Hence, the estimate of the entropy of counterion 
release given above should be considered a lower 
bound for the absolute entropy contribution. Other 
contributions not included in the theoretical anal¬ 
ysis apparently lead to experimental entropies that 
are higher by a factor of 2-3 (see cf. Table [^. The 
good agreement between the measured and the cal- 
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Table 2: The calculated standard binding free energy from the simulations AG^"” in comparison with the experimental 
ones for various salt concentrations and temperatures in units of kJ/mol. AG*"" is the direct output from the 
simulations which has to be corrected by AG‘^"" for the binding volume Vh to obtain the standard free energy of 
binding (see Methods) 


Conditions 

AG*™ (kJ/mol) 

Vt, inm^) 

AG^^"’" (kJ/mol) 

AG^™ (kJ/mol) 

^Qexp 

20 mM, 25°C 

-24.8 ±4.0 

3.0 

-1.5 

-26.3 ±4.0 

-27.9±0.2 

20 mM, 37°C 

-25.1 ±3.6 

5.9 

-3.1 

-28.2 ±3.6 

-29.2±0.2 

50 mM, 25°C 

-16.1±1.0 

16.0 

-5.7 

-21.8±1.0 

-23.5 ±0.5 

50 mM, 37°C 

-18.9±0.5 

7.5 

-3.9 

-22.8 ±0.5 

-26.2±0.2 


culated AG/,, however, demonstrates that these ad¬ 
ditional entropic contributions are canceled out by 
an enthalpic contribution of equal magnitude. This 
“enthalpy-entropy compensation“ is well-known 
for various processes as e.g. solute hydra tion, 

protein folding, or proteins association.^^Sl 

present comparison of theory and experiment al¬ 
lows us to discern among these terms and the lead¬ 
ing contribution to AG/,. 

Conclusions 

We presented a study of the binding of short PAA 
chains to HSA by combining calorimetry, that is, 
ITC and computer simulation. Both ITC experi¬ 
ments and simulation results show that there exists 
a strong attractive interaction between PAA and 
HSA at “the wrong side of pl“, where both are 
negatively charged. Computer simulations demon¬ 
strate that the binding of the PAA takes place at 
the Sudlow II site. ITC measurements for a se¬ 
ries of salt concentrations between 20 mM and 
100 mM show that the dependence of binding 
affinity AG/, on ionic strengths is mainly deter¬ 
mined by the counterion release mechanism and 
can be described by the Record-Lohman relation 
see eq. ?? (see Fig. |^. The binding affinity AG/, 
decreases with high ionic strength, until it practi¬ 
cally vanishes at around 750 mM. Both the anal¬ 
ysis of the experimental data as well as the simu¬ 
lations find that approximately 3 ions are released 
in the binding process. The binding affinity AG/, 
can be calculated from simulations with good ac¬ 
curacy. Combining simulations with calorimetry 
is hence a powerful tool to elucidate the interac¬ 
tion of proteins with given substrates or possible 
toxins. Thus, this combination of theory and ex¬ 


periment is now capable of solving biochemical 
problems of direct medical importance. 
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